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ABSTRACT 

In this paper, a controller order reduction method for linear parameter-varying systems is presented. 
The proposed method is based on the frequency-weighted balanced truncation technique, which 
has the advantage to reduce the order in a specific frequency range. The approach is discussed and is 
proved to preserve the closed-loop stability with a guaranteed upper error bound. Effectiveness and 
performance of the obtained reduced-order controller are investigated by applying it to an automo¬ 
tive semi-active suspension control. The obtained simulation results show that objectives such as the 
road handling and the passenger comfort realised with the reduced-order controller are kept in the 
same performance level as with the full-order controller. Moreover, a comparison with an other-order 
reduction method is shown and confirms the advantage of the developed method. 
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Notations 


R 

fields of real numbers 

-<(=<), ^(>=) 

negative (semi-negative), positive 
(semi-positive) definite 

Ko 

cone of n x n symmetric positive def¬ 
inite matrices 

n,o 

reachability and observability Grami¬ 
ans 

Pip), Qip) 

generalised reachability and observ¬ 
ability Gramians 

P n(p), Q n(p) 

frequency-limited generalised reacha¬ 
bility and observability Gramians 

Oi 

zth generalised Hankel singular value 

Yi 

zth frequency-limited generalised 
Hankel singular value 

S' 

complex conjugate transpose of 5 

A' 7 

(A- 1 ) 7 = (A 7 )" 1 

In 

n x n identity matrix 


1. Introduction 

Most of the physical systems are inherently nonlinear 
and almost all of them have parameter-dependent rep¬ 
resentations. Considering those aspects and the advan¬ 
tages of linear control theory, more and more control 
strategies use linear parameter-varying (LPV) models. In 
the same way, and as operating conditions may change, 


the closed-loop performance designed by an linear time- 
invariant (LTI) controller maybe degraded. To overcome 
this problem, the design of parameter-dependent con¬ 
troller is more suitable. Consequently, a major part of the 
recently developed control strategies is based on optimal 
and/or robust control. The -control strategies have 
the advantage to design controllers achieving stabilisation 
with guaranteed performance. However, these techniques 
usually produce high-order controllers. The design and 
the use of such high-order controllers can lead to numer¬ 
ical difficulties. This is why their use is still limited in the 
engineering field. Thus, for these practical reasons, low- 
order controllers are particularly preferred: their simple 
dynamics are easier to manage, they require less comput¬ 
ing cost and simpler software can be finally implemented. 
Then, an order reduction step appears of great interest. 
Simultaneously, the reduction process should always pre¬ 
serve the closed-loop stability and should guarantee a 
level of performance close to the one obtained with the 
full-order controller. In a general way, a low-order con¬ 
troller can be obtained through the direct method or 
the indirect ones (see Figure 1). Considering the direct 
method, the final low-order controller is obtained directly 
from the high-order plant and the order is generally fixed 
beforehand (Ankelhed, Helmersson, & Hansson, 2011; 
Apkarian & Gahinet, 1995). 

In the indirect ways, the low-order controller is 
the result of a model or a controller order reduction 
step: either a reduced-order plant is found and then 
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Figure 1. Controller order reduction scheme (Anderson & Liu, 1989). 


a controller is designed or a high-order controller is first 
designed and then a controller order reduction step is 
performed. The last method is widely used to obtain 
reduced-order controllers (Goddard, 1995; Kong, 2012). 
For LTI systems, a large investigation has been made on 
model order reduction procedures. In this domain, meth¬ 
ods based on balanced realisation are extensively used 
(Gugercin & Antoulas, 2004). First introduced by Mullis 
and Roberts (1976) and later in the systems and con¬ 
trol literature by Moore (1981), the balanced truncation 
(BT) has been a significant contribution to system the¬ 
ory. In Pernebo and Silverman (1982), the stability of the 
reduced-order model is preserved. Then, Glover (1984) 
has proposed an ^^-norm upper bound of the approxi¬ 
mation error. 

Based on this, the BT method has been proposed 
to reduce the order of controllers for LTI systems (Liu 
& Anderson, 1989; Zhou, D’Souza, & Cloutier, 1995). 
The extension of these reduction techniques to LPV sys¬ 
tems is still in progress. Generally, these methods substi¬ 
tute the use of LTI Gramians by using parameter-/time- 
varying equivalents (Sandberg & Rantzer, 2004; Wood, 
Goddard, & Glover, 1996). A generalised method with 
unbounded rate parameter model is given in El-Zobaidi 
and Jaimoukha (1998). In the same way, an effective BT 
method for J^-LPV controller order reduction is pro¬ 
posed in Widowati, Riyanto, and Saragih (2004). 

The BT produces good approximations in high 
frequency, whereas for control propose, an accept¬ 
able reduction is often required at intermediate and/or 
low frequencies. To solve this problem, the so-called 
frequency-weighted balanced truncation (FWBT) has 
been proposed first for LTI systems in Enns (1984). In 
this approach, two weighting functions are included in 
the procedure to reduce the model in a certain frequency 
range. However, the stability preservation is not guaran¬ 
teed anymore. To tackle this problem, a modified version 


assuming that the full-order model does not have any 
common pole/zero with the weighting functions is pro¬ 
posed in Lin and Chiu (1990). Then, Wang, Sreeram, and 
Liu (1999) and Varga and Anderson (2001) have devel¬ 
oped a new approach based on the Enns method without 
the previous assumption to guarantee the stability. It gives 
also an upper error bound. Other upper bounds have also 
been established by Sreeram and Anderson (1995), Kim, 
Anderson, and Madievski (1995) and Zhou (1995). The 
frequency range is explicitly defined instead of the two 
frequency weights in Gawronski and Juang (1990). More 
recently, in Gugercin and Antoulas (2004), the method is 
modified and an error upper bound to the relative error 
is given. Note that the stated techniques are developed 
for the stable LTI systems and no guarantee is given to 
preserve the passivity of the systems. Recent work in this 
direction are given in Li, Yin, and Gao (2014) and Li, 
Yu, and Gao (2015), where the ^^-norm of the error is 
bounded and the positive realness of the obtained model 
is guaranteed. 

Based on these recent results in the literature, the aim 
of this paper is to adapt the FWBT method (Gugercin & 
Antoulas, 2004) in order to reduce the order of an - 
LPV controller. For the development of this approach, 
the generalised Gramian framework is used (Shaker and 
Wisniewski, 2011). 

The paper is organised as follows: the LPV systems, 
their properties and the the full-order ^o-LPV con¬ 
troller design are introduced in Section 2. In Section 3. 
the generalised Gramian framework for LPV systems is 
given and the main proposed method is detailed and dis¬ 
cussed. To emphasise the effectiveness of the developed 
method, the obtained reduced-order controller is evalu¬ 
ated on a practical engineering problem: the control of 
an automotive semi-active suspension. The performances 
of the obtained reduced-order controller are investigated 
and compared with another reduced-order method based 
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on BT approach (Widowati, Riyanto, & Bambang, 2006) Proposition 2.1 (Briat, 2015): The system (1) is quadrati- 

in Section 4. A final conclusion is given in Section 5. cally stable if and only if there exists a matrix P e 0 such 

that the LMI 


2. Problem statement 
2.1 LPV systems 

LPV models represent systems whose state-space descrip¬ 
tions are known functions of varying parameters. This 
notion of LPV system has been introduced in the con¬ 
text of gain scheduling (Shamma & Athans, 1990) and has 
been also widely used in control design for nonlinear sys¬ 
tems (Packard & Kantner, 1996). 

Considering a compact subset A p C R s , an LPV sys¬ 
tem can be described by the following state-space realisa¬ 
tion: 


A T (p)P + PA(p) < 0 (5) 

holds for all p E A p . 

Moreover, the induced C 2 -gain notion is of a big inter¬ 
est for the proposed approach. Indeed, an upper bound is 
given to the C 2 -norm of error between the full-order and 
the reduced-order models. It can be defined as follows: 

Definition 2.2 (Induced C 2 -gain; Goddard, 1995): Given 
a quadratic-stable LPV system (1) with zero initial condi¬ 
tions, the induced £ 2 -gain is defined as 


x{t) = A(p(t))x(t) + B(p(t))u(t) > 

' y(t) = C(p(t))x(t) + D(p(t))u(t) 

( 1 ) 

where A : A p -> IT x ", B : A p -> M nxm , C : A p -> 
R pxn and D : A p -> R pxm . The state-space representa¬ 
tion (1) can be also written as 


j p 11*. 


i,2 = SUp SUp 


peV p ueC 2 II u II2 
u^0 


( 6 ) 


The induced £ 2 -norm represents the largest induced gain 
from inputs in jC 2 to outputs in C 2 over the set of all causal 
linear operators described by the LPV system. 


G{p) = 


'A(p) 

B(p) 

C(p) 

D(p). 


( 2 ) 


The exogenous parameter vector p is varying in a 
bounded set V p defined as 


2.2 LPV controller synthesis 

Consider a general LPV system represented under the fol¬ 
lowing state-space realisation: 


V p = {p : M + -> A p | pft) < pdf) < pft), i = 1 , 2 , ..., 5 } 
t pit) 


(3) 


x{t) = A(p(t))x(t) + Bi(p(t))w(t) + B 2 (p(t))u(t) 
z(t) = Ci(p(t))x(t) + Du (p(t))w(t) + D 12 (p(t))u(t) 'it 
y(t) = c 2 (p(t))x(t) + D 2 l (p(t))w(t) + D 22 (p(t))u(t) 

( 7 ) 


where Pi is the zth component of the p vector. 

The p function is not known in advance but is assumed 
to be bounded and measurable online. This assumption 
ensures that all trajectories of the system are contained in 
those ones of the LPV model. 

It is noted that some LPV models could be written as 
the function of system states (Sename, Gaspar, & Bokor, 
2013). Such models known as quasi-LPV systems are used 
for the modelling of nonlinear systems. 

Definition 2.1 (Quadratic stability): System (1) is said to 
be quadratically stable if the positive definite quadratic 
form 

V p :x\* x J P 0 x, P 0 € §: 0 (4) 

is a Lyapunov function for (1). A such Lyapunov func¬ 
tion is often referred to a common Lyapunov function or 
a parameter-independent Lyapunov function. 


where x(t) € R n , u(t) € R Hu , w(t) e R Uw , z(t ) € 

R Hz and y(t) € R n y are, respectively, the state, the 
input, the disturbance, the controlled output and the 
measured output. Then, A : A p —► R nxn , Bi : A p —► 
R nxn ^B 2 : A p -* R nxn \ Ci : A p -* R* xn , C 2 : A p -* 
D n : A p -* R nzXn -,D l2 : A p -> R HzXn \ D 21 : 

A p -* R n y xn ™ and D 22 : A p -> R n y xn \ 

The ,^ 0 -LPV controller associated to the LPV system 
(7) is defined by 

K( x x K (t) = A K (p(t))x K (t) + B K (p(t))y(t) 

{P) ‘ u{t) = C K (p(t))x K (t) + D K (p(t))y(t) 

( 8 ) 

where x K € R nK ,y(t) e and u(t) e are, respec¬ 
tively, the states, the inputs and outputs of the controller 
K(p). A k : Ap -> M nKXnK , B k : A p -> 1T kX % C K : 
Ap -> R n “ xnK and D K : A p -> R n » xn y. 
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The ,^ 0 -LPV controller synthesis concerns the 
design of an LPV global controller that guarantees 
both stability and performance for all parameter varia¬ 
tions defined in the set A p . To guarantee the closed- 
loop system quadratic stability and to satisfy the J4 ? 00 - 
performance criteria, the approach developed in Scherer, 
Gahinet, and Chilali (1997) is used to design the J^oo- 
LPV controller. It is assumed that 

• the matrices B 2 , D n , C 2 and D 2 \ are parameter inde¬ 
pendent and D 2 2 = 0; 

• the parameter-dependent matrix pair [A(p), B 2 (p)] 
is stabilisable and detectable VpG A p ; 

• the matrix [C 2 (p), D 2 fp)] has full-row rank Vp E 
A p . 

Then, for a given real positive scalar y and 
a parameter-dependent quadruplet matrices 
(A K , Bk, C k , Dr), there exist two parameter-independent 
symmetric matrices X and Y such that the sufficient con¬ 
dition that solves the ^o-LPV problem is given by the 
following LMIs Vp E A p : 


2.3 Full-order controller design 

Before introducing the approximation method, the full- 
order controller should be designed as shown in Sec¬ 
tion 2.2. Then, let 


T zw (p) = 


'A(p) 

B(p)" 

,C(p) 

D(p) 


= LFT 


'A(p) 

B(j>y 

C(p) 

Dip). 


'Ak(p) 


C K (p) 

D k (p). 


( 12 ) 


be the state-space realisation of the closed-loop sys¬ 
tem considering the augmented model (7) and its •K- 
controller defined in (8) with 


A(p) = 


A(p) + B 2 D K (p)C 2 B 2 C K (p) 
Bk(p)C 2 a k (p) 


B(p) = 


B\(p) + B 2 Dk(p)D 2 \ 

Bk(p)D 2 i 


C(p) = [Ci(p) + D l2 D K (p)C 2 D 12 C k (p)] and 
D(p) = D\\ip) + DuDk(p)D2i 


A(p)X + B 2 C K (p) + (*) 
AM+A\p)+ClDlip)Bl 
B\ip)+D\ x blip)Bl 
Q(p)X + D 12 Ck(p) 


* * 

TA(p) + Bk(p)Q + (*) * 

B|(p)7 + D^T(p) -yl n ^ 

Clip) + Di 2 byiip)C 2 Dnip) +Di 2 D k (p)C 2 

( 10 ) 


* 

* 

* 


-< 0 


(9) 


where terms denoted (*) are induced by symmetry, for 
example, 


~M + N + (*) *' 


~M + N + M J + N J K J ~ 

1 

_1 


-1 

_1 


The LPV controller of the form (8) is then reconstructed 
as 

D K (p) = £> K (p) 

(kip) = (C K (p) - D K (p)C 2 X)M- J 
Bk(p) = N~\B k (p ) - YB 2 D K (p)) 

Ak(p) = N~\A k (p) - YA(p)X - YB 2 D K (p)C 2 X)M T 
—B K (p)C 2 XM' J - N-'YBMp) (11) 

where M and N are defined such that MN J = I„ — XY 
and which can be solved through a singular value decom¬ 
position and a Cholesky factorisation. 


Note that the closed-loop system T zw ip) in Figure 2 is 
an LPV system in its own right which is quadratic stable 
by construction. In order to apply the proposed method, 
the T zw ip) state-space representation is assumed to be 
minimal. 

3. Main contribution 

In this section, the new controller order reduction 
method for LPV controllers is proposed. The method is 
based on the FWBT approach used for the LTI models 
order reduction. The FWBT is proposed to reduce the 
model order over a known frequency range [coi, co 2 ] in 
Gawronski and Juang (1990). This approach uses the clas¬ 
sical reachability and observability Gramians. 

Inspired by this work, the proposed method is an 
extension in two directions: 

• extended to the LPV case; 

• applied to reduce the controller order (rather than 
the model order). 
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Figure 2. Closed-loop scheme. 

To this end, the classical reachability and observabil¬ 
ity Gramians for LTI systems are introduced first. Then, 
the equivalent ones of the LPV case are defined. Based on 
this, the so-called generalised Gramians are introduced 
and then used to perform the proposed controller order 
reduction approach. The relation between the Gramians 
and generalised Gramians is given in Lemma 3.1. 

3.1 Frequency-limited generalised gramians 

The proposed method is based on the FWBT method 
where the key notions are the so-called reachability and 
observability Gramians. 

Definition 3.1 (Gramians): Considering the quadratic 
stable and minimal realisation of the closed-loop LPV 
system T zw (p) given in (12), Let us define 7 Z(p) and 
O(p), the related reachability and observability Grami¬ 
ans, respectively; they are defined as the solution of the 
two parameter-dependent Lyapunov equations 

- V ». + A(p)TZ(p) + U(p)A J (p) 

tr 9 a- 

+B( / o)B T (p) = 0 (13) 

y-v + x J (p)0(p) + (D(p)A(p) 

tr d p> 

+ C T (p)C(p) = 0 (14) 

For a given parameter trajectory p , let O p (t, 0) be 
the state transition matrix. Then, the functional IZ(p) 
and 0(p ), respectively, the reachability and observabil¬ 
ity Gramians of the closed-loop LPV-system (12), are 
expressed such that 

r° 

nip) = / <E> P (0, OB(p(O)B T (p(O)^(0, t)dt 

J—o o 

(15) 


O(P) = / <jT(f, 0)C T (p(0)C(p(0)4> P (f, o )dt 
Jo 

( 16 ) 

Definition 3.2 (Generalised Gramians): Let P(p) and 
Q(p) be, respectively, the generalised reachability and 
observability Gramians of the quadratic stable and mini¬ 
mal closed-loop LPV system (12) given as 
V(p, v) G A p x A v 

A dP (p) t 

- V °i -1- A(p)P(p) + P(p)A t (p) 

tr 9 ^- 

+ B(p)B T (p) -< 0 (17) 

V >p d( ^ (p) + A T (p)Q(p) + Q(p)A(p) 

G 3 Pi 

+ C T (p)C(p) -< 0 (18) 

Indeed, Lyapunov inequalities (rather than Lyapunov 
equations) are solved to compute generalised Grami¬ 
ans. This linear matrix inequality (LMI) approach to 
the model reduction problem is particularly useful when 
some structures need to be preserved in the process of 
model reduction. Controller reduction is a typical exam¬ 
ple of this type of problems (Shaker & Wisniewski, 2011). 
Note that the physical interpretations of generalised 
Gramians are similar to ordinary Gramians. Considering 
Tzw (p)i'R'(p) and O(p), the following lemma introduces 
useful results about the relation between Gramians and 
generalised Gramians. 

Lemma 3.1: Let T zw (p) be a minimal state-space realisa¬ 
tion of the quadratic stable and minimal closed-loop system 
defined by (12). Then, Vp € V p 

H(p) < P(p(0)) 

O(p) < Q(p(0)) 1 ; 

where TZ(p), (9(p), P(p) and Q(p) are the solutions of 
(15), (16), (17) and (18), respectively. 

Proof: By multiplying on the left by 0^(0, t) and on 
the right by 0^(0 ,t), the relation (17) yields 

s gp 

- O p (0, 0y>,— 4> p (0, t) + <D P (0, f)A(p(0)P(p)<J>^(0, t) 

tr d pi 

+ $ p (0,f)P(p)A T (p(0)^(0, t) 

+ 4> p (0, f)B(p(f))B T (p(f))4> p (0, t ) .< 0 


( 20 ) 
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Definition 3.3 (Frequency-limited Gramians): Let 
7£^(p) and Oq(p), two frequency-dependent terms, 
defined as 


Using the fact that t 0 ) = -<F p (t, t 0 )A(p(t 0 )), 

then (20) yields 

-^(4> p (0,OP(p)<l>p(0,O) 

+ 4>p(0, t)B(p(t))B J (p(t))^ J p (0, 0 -< 0 (21) 

which on integrating over the semi-infinite time axis 
(— oo, 0] and considering that lim t —— 00 4> p (0, f) = 0, 
gives 

P(p(0)) >- / <t> p (0,t)B(p(t))B J (p(t))<t> J p (0,t)dt 

J —00 

K(p) 

( 22 ) 

Similarly, let us multiply on the right by the closed-loop 
state transition matrix <$> p (t, 0) and on the left by (t, 0), 
the relation (18) 

0 } ±A p (t, 0) + 0 T (t, 0)A T (p(t))Q(p)<P fi (t, 0) 
+<t> T p (t, 0)Q(p)A(p(f))<t> p (f, 0) 

+<t> J p (t, 0)C T (p(f))C(p(f))<I>p(f, 0) -< 0 (23) 

Using the fact that ^<t> p (t, to) = A(p(t))<F p (t, to), then 
(23) yields 

j t (4> T p (t,0)Q(p)<t> p (t,0)) 
+<t> T p (t,0)C T (p(t))C(p(t))<t> p (t,0) ^0 (24) 

Integrating the last expression from 0 to +oo gives 

r+oo 

Q(p(0))> / <t>l(t,0)C T (p(t))C(p(t))<i> p (t,0)dt 

Jo _^ 

O(p) 

(25) 


3.2 Frequency-limited order reduction 

For a given parameter trajectory p , let consider O p (f, 0) 
the state transition matrix of the closed-loop system (12). 
Then, let us define 

fp(t) = O p (0, t)B(p(t))H(-t) 

g p (t)= C(p(t))O p (t,0)H(t) K } 

where H is the Heaviside step function. 

By considering £2 = [coi;co 2 \ the frequency range 
where order reduction is desired to be better, the follow¬ 
ing definition is given. 


Tln(p) = T^co 2 (p) ~ Kco! ( P ) and Oq ( p) 

= 0*(p)-0» 1 (p) 


where 


/ +CD 

F p (w)F;(w)dw (27) 

-CD 

/ +(D 

G*(w)G p (w)dw (28) 

-CO 

with F p (respectively G p ) is the Fourier transform of 
f p (respectively g p ). Then, the functional 7 Zq(p) and 
Oq(p), respectively, the frequency-limited reachability 
and observability Gramians of the closed-loop LPV sys¬ 
tem (12), are defined as the solutions of 
V(p, v) E A p x A v 

- + A(p)Hn(p) + Hn(p)A J (p) 

tr 9 ^' 

+ n Q (p)= 0 (29) 

y-v aOft(p) + A J (p)dn(p) + On(p)Mp) 

tr 9 P' 

+O q (p) = 0 (30) 

The quantities 7 Z&(p) and Oq(p) have the following 
eigenvalues decomposition: 


ft-n(p) = Uj 2 (p)diag(A 1 (p),A.„ K (p)) l£(p) 

(31) 

Oq(p) = V Q (p)diag (Siip),..., S„ K (p)) V^(p) 

(32) 

with |Ai| > • • • |A„ k | > 0 and \8i\ > • • • |S„ K | > 0. Let u p 
< riK and v p < n k be, respectively, the ranks of 7 Zq(p) 
and Oq(p). Based on these definitions, let us define the 
two quantities 


Bfi(p) = l/ n (/o)diag(|Ai(p)|5,..., \X Up (p)| 0,..., 0), and 

(33) 

Cfi(p) = diag(|5i(p)|U ., |5„ p (p)|5,0,..., O)V^(p) 

(34) 







INTERNATIONAL JOURNAL OF CONTROL @ 2037 


Definition 3.4 (Modified frequency-limited generalised 
Gramians): Consider Kip), the full-order stabilising LPV 
controller given in (8). Let P^(p) = p ( p) j Oand 

Qciip) = ^ 1 Q (p) q° ( p) j ^ 0 be the modified frequency- 
limited reachability and observability Gramians defined 
as the solutions of the following Lyapunov equations 
V(p, v) E A p x A v : 

- J2 + A(p)P Q (p) + P Q (p)A T (p) 

i^i dpi 

+ Bq(p)Bq(p) = 0 (35) 

J2 + A T (p)Qo(p) + Q n (p)A(p) 

tr dpi 

+ = 0- (36) 

For the generalisation, we have the following inequalities: 


- T2 v i —f — + A(p)Pn(p) + Pj2(p)A T (p) 

tr dpi 

+ B^(p)BT(p) ^0 (37) 


with^i(p) > § 2 (p) ^ ... > ?„(p),and 

T 2 - 1 (p)P 2 (p)T 2 - t (p) 

= T 2 t (p)Q 2 (p)T 2 (p) 

= s 2 (p) 

= diag(yi(p),..., y r (p), Yr+\(p), ■ • •, Vn K (p)) 

^21 ^22 

(40) 


with Ki(p) > y 2 (p) > • • • > Kr(p) > Kr+i (p) > 

y r+2 (p) > • • • > Kn K (P) are the frequency-limited 
generalised Hankel singular values of K(p) and r is the 
desired order for the reduced-order controller. 

The balanced realisation of K(p) can be written as 


K(p) = 


T 2 - 1 (p)A K (p)T 2 (p) 

T 2 - 1 (p)Bk(p)' 

Ck(p)T 2 (p) 

Dk(p) 


^k(p) 

Bk(p)' 

Qc(p) 

Ac(p) 


(41) 


Further, K(p ) is partitioned as conformably with £ 2 (p) 
as 


V + A T (p)Qn(p) + Q n (p)A(p) 

tr 9pi 

K 0. (38) 




^k(p) ^l K i 2 (p) 

r ? Kip) " 

K(p) = 


^K2i(p) ^K22 (p) 

BK2(P) . 



[Qc(p) Ck2(p) 


DM 


with P Q (p) = P 2 ( P) 

[ Qi(p) o 1 o 

L 0 q 2 (p)J > u - 


>~ 0 


and Q^(p) = 


If the block diagonal solutions P^(p) and Q^(p) exist, 
then let Tfp) and T 2 (p) be two non-singular matrices 
given such that 


T 1 - 1 (p)P 1 (p)r i - T (p)=T 1 r (p)Q 1 (p)T 1 (p) 

= S,(P) (39) 

= diag(§i(p), ...,§„(p)), 


Finally, a truncation step is performed to obtain a 
reduced-order controller. 

Definition 3.5: Given the balanced realisation K(p) 
defined in (42), let K(p) be the truncated realisation to 
the rth order and denoted as follows: 


K(p) = 


'am 

Bk(p)' 

Qc(p) 

O k (p) 


(43) 


Furthermore, the reduced-order parameter-dependent 
closed-loop system is given as 


Tzw(p) = 


A(p) 

B(p)' 

C(p) 

D(p) 


A(p) + B 2 D K (p)C 2 5 2 C K (p) 
£k(p)C 2 A k (p) 


Ci(p) +T ) i2^ ) k(p)C 2 D J2 Ck(p) 


Blip) + B 2 D K (p)D K (p)D 2 i 
Bk(p)D 2 i 


Du ip) + T)i 2 Dk ip)D 2 \ 


( 44 ) 
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Theorem 3.1: Suppose K(p) is the stabilising parameter- 
dependent controller defined in (8) such that the closed- 
loop transfer T zw (p) defined in (12) is minimal, quadratic 
stable, and there exist Lyapunov inequality solutions Pq(p) 
and Qq(p) such that (37) and (38) are satisfied. Let 
K(p) be the reduced-order controller defined in (43) and 
obtained by truncation. Then, the closed-loop system with 
the reduced-order controller T zw (p ) defined in (44) is sta¬ 
ble. If, in addition, 

rank[B(p), B Q (p)] = rank[B n (p)] (45) 

and 

rank[C 7 (p), C£(p)] = rank[C£(p)] (46) 

then, f zw (p) is quadratic stable and satisfies 

n K 


Since 7 Z^(p) and 0&(p) are not guaranteed to be 
positive definite, the stability of the reduced-order con¬ 
troller is not guaranteed. Then, an idea based on eigenval¬ 
ues decomposition presented in Gugercin and Antoulas 
(2004) is proposed here to guarantee the stability by pro¬ 
viding an upper error bound. Indeed, the solution of (37) 
and (38) is performed instead of (29) and (30). In addi¬ 
tion, the modified frequency-limited Gramians in (35) 
and (36) are used instead of the ordinary ones defined in 
(35) and (36). In fact, Lemma 3.1 shows that the gener¬ 
alised Gramians could be an approximation of the ordi¬ 
nary Gramians. The transition to generalised Gramian 
framework might induce less accurate approximation 
but the order-reduction error still bounded. Finally, the 
expression of the upper error bound given in (47) is found 
according to the following. 

Let T zw (p ) be the frequency-limited balanced realisa¬ 
tion of the full-order closed-loop system. Then, T zw ( p ) is 
defined by 


II T ZW (P) 


T zw (p)\\i,2 < 


2||Mp)IIJ/c(p)IL 


E 

i=r +1 



T — 

r T~ l (p)A(p)T(p) 

T- 1 (p)B(p) 

(47) 

J-zw \P ) — 

[ C (p)T(p) 

D(p) 


where J B (p) := diag(|Ai | - ^ (p),1^1“^ (p), 0,..., 0) 

T A(P) A U (P) 

i 

r B(p) ■ 

Uq(p)B(p) and = 

[a 2 i(p) A 22 (p) 

i 

|b 2 (p) 

7c(p) := C(p)Vj 2 (p)diag(|<$i| - ^(p),..| K p \~Hp), 

r - - i 



o,.. •, 0). 

[C(p)C 2 (p) 


D(p) 


Proof: The reachability and the observability Grami- 


ans given in (15) and (16) can be expressed as 

p+OG 

n(p) = / fp(j)f*(j)dr 

J —00 

(48) 

poo 

O(p) = / g*(r)gp(r)dr 

J —00 

(49) 

where/ p (r) andg p (r) are given in (26). 

Then, using Parseval relationship, the reachability and 

the observability Gramians could be expressed as 

follows: 

1 f +0 ° 

K(p) = — F p (w)F*(w)dw 

27r J — qq 

(50) 

1 f +0 ° 

O(p) = — / G (w)G p (w)dw 

J _ 00 

(51) 


By considering F p (resp.g p ), the Fourier transform of 
f p (respectively^), we can define 7 Z M (p) and C^(p) as 
the limited reachability and observability Gramians given 
in (27) and (28). Then, by setting Q = [coi ;co 2 ] the fre¬ 
quency range where approximation is desired to be bet¬ 
ter. Then, new modified terms TZ^(p) and Oq(p) are 
defined. 


where T(p) = diag(Ti(p), T 2 {p)). 

By considering the assumptions (45) and (46), there 
exist J B (p) and J c (p) such that B(p) = B n(p)J B (p) and 
C(p) = 7c(p)C^(p). On the other hand, the reduction 
error is expressed as 

II T ZW (P) - tw(p)\\i ,2 = ||C(p)(s7- A(p))“ 1 B(p) 
-C(p)(s7-A(p))- 1 B(p)|| / , 2 
= II/c(p)[Cq(p)(s7 — A(p)) -1 Bq(p) 

-C(p)(57-A(p))- 1 B(p)]./ B (p)|| i , 2 

n K 

<2||/b(p)||J|/c(p)IL J2 

i=r +1 


Remark 3.1: Assumptions (45) and (46) mean that there 
exist Jb(p) and J c (p) such that B(p) = B n(p)J B (p) and 
C(p) = Jc(p)Cq(p). In addition, by following the steps 
in Wang et al. (1999), Anderson and Clements (1981) and 
Imran, Ghafoor, and Sreeram (2014), it was shown that 
Assumptions (45) and (46) are almost always true. Hence 
we expect that our approach will apply in most of the 
cases. Indeed, during our simulations, the assumptions 
have always been satisfied. 
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Algorithm 1 J^-LPV controller order reduction 
Considering the LPV plant described in (7), the reduced- 
order controller can be computed as follows: 

Inputs: (A(p), B(p), C(p), D(p)). 

Outputs: (A(p), B(p), C(p), D(p)). 

Assumptions: (A(p), B(p), C(p), D(p)) minimal. 
Algorithm: 

(1) Compute K(p) the full-order controller according 
to Section 2.2. 

(2) Compute the closed-loop system T zw given in (12). 

(3) Compute the generalised Gramian P(p) and Q(p) 
solutions of (17) and (18), respectively. 

(4) Compute the balanced realisation K(p) of the full- 
order controller K(p) by the following: 

(a) Find T 2 , the basis change matrix according to 
(40) (for instance, Procedure in []Moorel988 
can be used). 

(b) Compute the balanced realisation K(p) as 
defined in (42). 

(5) Compute the reduced-order controller K(p) from 
K(p) by truncation. 


4. Case study: semi-active suspension control 
4.1 LPV model 

In this section, the given method to reduce the LPV con¬ 
troller is applied on a semi-active automotive suspen¬ 
sion presented in Do, Sename, Dugard, and Boussaad 
(2011). Actually, when suspension modelling and con¬ 
trol are considered, the vertical quarter car model is often 
used. This model allows to study the vertical behaviour of 



a vehicle according to the suspension characteristic (pas¬ 
sive or controlled). Figure 3 shows the so-called vertical 
quarter car. Then, the dynamical equations of the system 
are given by 

m s Z s — ^s^def Fmr (53) 

^us^us — ^s^def + Fmr ^tGus ^r) 


where F mr is the magneto-rheological force generated by 
the semi-active suspension. According to the nonlinear 
model of Guo, Yang, and Pan (2006), F mr can be expressed 
as follows: 


F m r — ^2 ^Zdef T ^ £def^ T tanh ^#3 ^Zdef T ^ £def^^ 


(54) 


with Zdef = z s — z us is the damper deflection (must be 
measured or estimated) and z<± e f = — z us is the damper 

velocity. Parameters a 2 , a 3 , v 0 and x 0 are constant, and a\ 
is the controllable force such that a x E [ai min ; ai max ]. 

By defining 

Pi = tanh (a 3 (z de f + ^z de f)), 
c mr = a 2 : MR damping coefficient, 
k mY = a 2 — : MR stiffness coefficient, 

Then, a state-space representation can be given by con¬ 
sidering the state vector x s = [Zs^Zus^us] 1 and the exoge¬ 
nous input w = z r , as follows: 


*s — A s x s + B s pid\ + B S \W 
• £ = C s iX s + D s ipiai (55) 

k y = c s x s 

with 


/ 0 


Ac = 


B s = 


^s+/c mr C mr 

ks “T k mr 

m s m s 

m s 

0 

0 

0 

l k s +k 
' m u; 

mr C mr 

k s +k mr - 

s m us 

m us 

0 ^ 


0 \ 

-1 

5 B S 1 = 

0 

m s 

0 

0 

\£) 


\ m us / 


° \ 

Cmr 

m s 

1 

_ Cmr / 

m us ' 


The measurement output is y = z s — z us , and the con¬ 
trolled outputs are chosen as z = [z s z s ] J , respectively, 
the acceleration and the displacement of the sprung mass. 
Then 


Gi 


_ ks ~F^mr C m r ~F^mr C m r \ 

m s m s m s m s | 

1 0 0 0 / 


Dsi 


m s 

0 


,c s = (10-10). 


Figure 3. Quarter-vehicle model. 
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However, two constraints must be satisfied: 

(1) The control signal a\ must be positive (dissipative 
constraint). 

(2) The input matrices B s p i and D s ipi must be con¬ 
stant to satisfy the LPV-^^o synthesis assump¬ 
tion. 


The passivity problem is solved by defining a new con¬ 
trol signal u — a\ — T 0 , where T 0 is the mean value of a\ 
(T 0 = (ai max — 0imin)/2). Then, the problem of the pas¬ 
sivity on ai is recast to a simple saturation problem on u 
(wE [ — Tq ; Tq]). With these modifications, (56) yields 

. Pi 

X s = (As + B s 2 ~ -C S 2)x s + B s p\U + B S \W 

G s 2*s 

z = C sl x s + DsiPiU ( 56 ) 

y — C s x s 

where 

Bs 2 = (o -£ 0 ^) T andC s2 = «3 -*3? 

To overcome the second problem, Apkarian and Gahinet 
(1995) propose to add a strictly proper filter T to 
make the controlled output matrices independent of the 
scheduling parameters 


T : 



A f Bf 
C f 0 



(57) 


Then, by defining p 2 = and v = (x s Xf ) T , the sys¬ 
tem (56) can be represented as 

x — A(pi, Pi)x + Bu c T B\W 
• z = Ci(p u p 2 )x (58) 

y — Cx 


where 


A(p u Pi) = 
B = 


(A s + PiB s iC S 2 P\B s Cf 



Ci(pi, Pi) = (Csi PiD sl Cf) and C = (C s 0) 


4.2 J^-controller synthesis 

By considering the LPV model (58), an J4 ? 00 -controller 
is designed to guarantee the internal closed-loop stabil¬ 
ity and to satisfy some required performance. In fact, the 
main objective and challenge of a controlled suspension 


system is to improve the comfort for car passengers simul¬ 
taneously to the performance on road holding. The pas¬ 
senger comfort can be improved by isolating the vibra¬ 
tions transmitted from the road surface. Then, the ‘fre¬ 
quency response from the road profile z r to the vehi¬ 
cle chassis acceleration z s must be kept small in the low- 
frequency range. Then, a weighting function is designed 
as 


JAr s 2 + §n wns + w\ x 

Wz = Wz - - ~y~ 

S s 2 + 2^\ 2 ID\ 2 S + w l2 

Furthermore, the road holding is evaluated from the 
unsprung mass (wheel) oscillations with respect to the 
road profile. This transfer should be kept small at high 
frequencies. Then, V\4 us is designed as 

JAT S 2 + § 21 ^ 21 $ + w>2\ 

W z = w z — - 

US S 2 + 2^21^>2lS + ID 22 

W Zr = 5 x 10 3 is the road profile gain. Finally, the filter 
introduced in (57) is given as T = it i s designed 
with a large bandwidth to decouple the input and the 
varying parameters, where 


Wz s = 1, §n = 0.1, §12 = 1, inn = 2 tt x 1 rad.s \ 
w\ 2 = 27r x 3 rads -1 

w Zus = 10, §21 = 0.3, §22 — 1, id 2 \ = 2 tt x 9rads _1 
id22 = 2 tt x 9 rads -1 , Wf = 90.34 


Then, an interconnection between the LPV model and 
these weighting functions are presented in Figure 4. 

Model parameters are obtained by considering experi¬ 
mental data in Tudon-Martinez et al. (2013) and given as 
Table 1. 



Figure 4. ^ > OQ —control scheme. 
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Table 1. Values of (rank[B(p), B^(p)]/rank[B^(p)]; rank[C T (p), Cj 2 (p)]/rank[Cj 2 (p)]) for the several 
frozen values of [p h p 2 ). 


p,/p 2 

- 1.0 

-0.5 

0.0 

0.5 

1.0 

0 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14 ; 14 / 14 ) 

0.25 

(14/14; 14/14) 

(13/13; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14 ; 14 / 14 ) 

0.5 

(11/11; 14/14) 

(13/13; 14/14) 

(11/11; 14/14) 

(12/12; 14/14) 

(11/11; 14/14) 

0.75 

(12/12; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14 ; 14 / 14 ) 

1.0 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14; 14/14) 

(14/14 ; 14 / 14 ) 


To carry out a controller satisfying these objectives, 
the ^^o-LPV synthesis is designed by using solution for 
polytopic systems: it consists in finding a global LPV con¬ 
troller K(p i, p 2 ) which is a convex combination of local 
controllers obtained by solving the LMIs set at each ver¬ 
tex (formed by limits values of the varying parameters). 
All varying parameters are bounded: p\ E [ — 1; 1], p 2 
E [0; 1]. 

For more details and explanation on ^^o-LPV control 
synthesis, see Apkarian and Gahinet (1995) and Scherer 
et al. (1997). The design method for LPV systems is used 
like in Apkarian, Gahinet, and Becker (1994). 


43 Numerical issue 

The proposed method requests the solution of two Lya¬ 
punov inequalities with an infinite number of constraints. 
These sets of infinite LMIs can be solved by gridding 
techniques. Then, some approximations must be made by 
gridding the set A p with finite number L of points {p z }[ =1 
(Lee, 1997). Moreover, the infinite variables P^(p) and 
Q^(p) in LMIs (37) and (38) are approximated by com¬ 
binations of scalar basis functions such as 

N P 

p q(p) = X>iG°) p ; v 0 

7=1 

Nq 

Qq(p) = X>(p)Q, >- 0 

j=l 

where Pj = Pj, Qj = Qj. There is a large freedom in the 
choice of basis functions (Wood, 1995). For this example, 


the following choice is made: 

1 = {<Pj})L 1 = { 1 , Pi, Pi, P\, pj, PlPl, P\P\, p\pl, 
P\P\, P\, P\, P\P\, Pi Pi } 

The main consequence of this approximation is that the 
number of LMIs to be solved is finite and is 2L(2 S + 1 + 1) 
where s is the number of parameters, i.e. s = 2. 

The full-order controller is designed using the proce¬ 
dure developed in Scherer et al. (1997). Then, an eighth- 
order controller Kf u \\ is obtained. The proposed method 
is used to reduce Kf u \\. Then, Kfwbt is produced. 

To test its effectiveness, a comparison with the method 
developed in Widowati et al. (2004) (LPV BT), namely 
Kbt, is performed. The feasibility problems (37) and (38) 
are convex. Using Matlab LMI Control Toolbox, con¬ 
trollers are reduced to the fifth-order obtained heuristi- 
cally by trial-and-error approach. Therefore, a frequency 
and time analysis are performed. 

5. Results and discussion 

The first evaluation is represented in Table 2 where first 
the assumptions (45) and (46) are checked. Note that for 
25 frozen values of (pi, p 2 ), the rank assumption is sat¬ 
isfied for all these points. This fact confirms Remark 3.1 
and allows us compute an upper bound. Precisely, Table 3 
expresses the upper bounds and norm of the error 


Table 3. Parameter values. 


Parameter 


Value 

Sprung mass 

(m s ) 

470 

[kg] 

Unsprung mass 


110 

[kg] 

Tyre stiffness coefficient 

(k t ) 

86378 

[N/m] 

Spring stiffness coefficient 

oy 

270000 

[N/m] 


Table 2. Values of the || T zw (p) - T zw (p) || ij2 and the upper bound for several frozen values of (p h p 2 ). 


pM 

- 1.0 

-0.5 

0.0 

0.5 

1.0 

0.0 

1.36 < 5.07 

1.35 < 3.32 

1.33 < 4.08 

1.31 < 3.37 

1.30 < 3.72 

0.25 

0.32 < 9.88 

0.32 < 9.93 

0.32 < 10.38 

0.31 <10.39 

0.31 < 9.80 

0.5 

0.00 < 1.79 

0.00 < 0.13 

0.00 < 2.22 

0.00 <11.81 

0.00 < 2.03 

0.75 

0.31 < 9.97 

0.31 < 9.92 

0.30 <10.58 

0.30 <9.86 

0.29 < 9.89 

1.5 

1.22 < 3.77 

1.20 < 3.04 

1.19 < 4.46 

1.18 < 2.88 

1.17 <2.90 
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Z s jZ r 



Figure 5. Bode diagram magnitude of the transfer T z . 

s r 


Z s /Z r 



Figure 6. Bode diagram magnitude of the transfer T ±z . 

s r 
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^US /' Z r 



Frequency (Hz) 


Figure 7. Bode diagram magnitude of the transfer T z z . 

us r 


values. These results evince that the upper bound is cor¬ 
rectly positioned (the gap is positive). However, we note 
that this upper bound is not tight to error Hoo -norm. 

5.7 Frequency analysis 

The Bode diagrams at several frozen values of pi and P 2 
(25 points) of the three transfer functions T Zs , T Zs and T Zus 
are shown, respectively, in Figures 5-7. In fact, the fre¬ 
quency behaviours of the chassis position z s and accel¬ 
eration z Y are chosen to be analysed in order to observe 
the comfort performance regarding the road profile input 
z r . The wheel position signal z us is also analysed to test 
the road holding. Then, the weighting functions W Zs and 
W Zs designed in Section 4.2 limit the amplification of the 
previously cited transfers in low-frequency range (around 
[1 ; 10] Hz). In fact, the human sensitivity to verti¬ 
cal vibrations is important in this frequency range (Do 
et al., 2010). For this reason, the frequency interval of the 
proposed frequency-limited FWBT method is chosen as 
[1; 8] Hz. 

In Figure 5, note the reduced-order closed-loop sys¬ 
tem produced by FWBT approximates well the full-order 
closed-loop system in the chosen frequency range [1 ; 
10] Hz. In this interval, the reduced-order closed-loop 


system produced by the unweighted BT fails. In fact, an 
important gap appears as 2 and 3 Hz which is exacted 
as the BT is known to guarantee good approximation at 
high frequency. The same comment is given in Figure 6 
where FWBT fits the full-order closed-loop system in all 
the shown range unlike BT method that misses the peak 
around 2 (a resonance frequency). These results are more 
explicit when observing Figure 7. Indeed, FWBT gives 
a good approximation when BT fails (2 and 8 Hz). The 
other important fact stated by these results is the sensitiv¬ 
ity against the parameters variation. Actually, a dispersion 
of pi and p 2 values is induced by the reduction step. This 
degradation is expected since the given application is a 
qLPV system where p x and p 2 are depending on the state 
vector. Then, every decrease in the order (states) affects 
the pi and p 2 values. However, this loss when reducing 
is under control for two reasons: the first one, the sta¬ 
bility of reduced-order closed loop is preserved and the 
error is guaranteed limited. The second one is that this 
dispersion is weak in the required frequency range. Thus, 
performance is not affected. 

Remark 5.1: As Bode diagram for nonlinear systems is 
not possible, a pseudo-Bode plot is then proposed in 
Poussot-Vassal (2008) 
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Chassis position 



Time [s] 


Figure 8. Time response on a bump of 0.01 m x2 m. 


5.2 Time analysis 

In the time domain, the several controlled suspensions 
are travelling a bump of 0.01 m x2 m for a vehicle speed 

8.3 m/s (i.e. 30 km/h). 

It is observed that the time response confirms the 
contribution of the FWBT reduction method. In fact, in 
Figure 8, the chassis is stabilised rapidly (1 s after the 
perturbation) without overtake on the suspension unlike 
the suspension with K^j. This observation preserves the 
required comfort’ performance. Moreover, for z us , the 
wheel equipped with a KWbt suspension keeps almost 
the same profile of the road although its variations which 
which respects the road-handling performance. The K^t 
suspension generates an infinite perturbation just after 
the bump (after 1 s). Note also that temporal test draws 
two output signals (chassis and wheel positions) regard¬ 
ing the input (the road profile) and by the way there will 
be just one plot of each transfer besides the several plots 
in the frequency responses. 


6. Conclusion 

In this paper, the problem of LPV controller order reduc¬ 
tion is investigated, and the approximation in a limited 
frequency range is considered, based on a similar work 
on the LTI case, a new frequency-limited method for 
LPV controller order reduction is derived. The obtained 
reduced-order controller is proven stable and the degra¬ 
dation in the closed-loop performance is guaranteed 


bounded. An application of the proposed method to 
reduce the order of a semi-active suspension controller 
is performed. The obtained results show its effective¬ 
ness to preserve the closed-loop stability and to respect 
the required performance specifications. The compari¬ 
son with another existent LPV order reduction technique 
based on the the classical BT, confirms the contribution 
of the proposed method. 
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